Dynamic scaling theory of the Kosterlitz-Thouless-Berezinskii transition: ubiquitous 

finite size effects 
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A numerical study of the neutral two-dimensional (2D) lattice Coulomb gas is performed to examine 
dynamic scaling, finite size effects, and the dynamic critical exponent z of the Kosterlitz-Thouless- 
Berezinskii transition. By studying large system sizes (L = 100), we show z — 2.0 ± 0.2 using 
Fisher-Fisher-Huse (FFH) scaling. We also present evidence that the vortex correlation length is 
finite below the transition temperature, in contrast to conventional wisdom. Finally, we conclude 
that previous findings for a variety of experimental systems that z — 5.6 using FFH scaling indicates 
that finite size effects are ubiquitous in 2D superconductors and Josephson Junction Arrays. 



I. INTRODUCTION 

The universality class of the two-dimensional (2D) XY 
model has been fascinating from the start with its phase 
transition being driven by topological excitations, the un- 
usual, exponential temperature dependence of the corre- 
lation length, and its apparent presence in systems rang- 
ing from 2D superfluids and superconductors to liquid 
crystals and magnetic systems.ErO Despite its unique- 
ness, the detection-p£the Kosterlitz-Thouless-Berezinskii 
(KTB) transitionJlnj characterized by an unbinding of 
vortex pairs, has been challenging, but several techniques 
have been developed to probe its critical behavior. In 
two-dimensional superconductors and Josephson Junc- 
tion arrays (JJA's), an examination of the dc transport 
properties in zero magnetic-field has been the primary 
method, but the magneto-resistanceu and the ac kinetic 
inductanceB have also been used. More recently, mag- 
netic flux noise measurementsQ have been used to study 
the KTB transition. And, as we discuss shortly, numeri- 
cal methods and scaling approaches have also been used. 
The advantage of the scaling methods is that they can 
reveal information on the dynamic critical exponent z, 
one of the key exponents in the KTB critical behavior, 
which relates the relaxation time r to the vortex corre- 
lation length £ via r oc £ z . 

While there are many methods for studying the KTB 
transition,!] the most common is the "conventional" 
method which involves an analysis of the dc voltage- 
current (V-I) isotherms and the temperature dependent 
resistance R(T). In this approach, the transition temper- 
ature T c is determined via the V-I exponent a(T, I) (V oc 
I a ), which is predicted to have a value of 3 at T c (assum- 
ing z = 2,) and whose value should decrease rapidly for 
weak currents near T cn The Ohmic resistance is then 
fitted by the Kosterlitzcl or Minnhagen formEl for the re- 
sistance, which are respectively, R oc exp[— bzj \JT — T c ] 
and R oc exp[-6z^(T c0 - T C )/(T - T c )], where b is a 



system-specific parameter and T c q is the mean-field tran- 
sition temperature. 

There are many weaknesses to the conventional 
approach.!] For example, it is valid only in the weak- 
current limit where few experiments can probe effectively 
and where finite-size effects dominate. Moreover, the ex- 
ponent a has a current dependence and so determining 
its value uniquely for a given temperature is not possi- 
ble. In fact, researchers can use this flexibility to their 
advantage to show a rapid variation of a with tempera- 
ture near T c . This is done by determining a in different 
current ranges for each temperature so as to achieve the 
desired behavior. (The most accurate method would be 
to determine a for the V-I isotherms for the, same cur- 
rent so as to represent a single length scale.B) Another 
problem with the traditional approach is the number of 
fitting parameters in the resistance formulas, which al- 
low excessive freedom for interpreting results. For these 
reasons, it is important to examine other approaches to 
verifying and studying KTB behavior. _ 

The Fisher-Fisher-Huse (FFH) dynamic scaling,!!! 
which says that V£ z /I is a function of I£,/T, has its 
advantages over the conventional approach. The most 
prominent is that it is valid over an extended current 
range and therefore incorporates the current dependences 
of a(T, I). Moreover, FFH scaling can be used to deter- 
mine z. Recently, it has been used to study experimental, 
dc V-I data but yields a large value of the dynamic crit- 
ical exponent .z ~ 5.6 for a wide variety of samples and 



systems. I 



It has recently been reiterated that the 



FFH scaling is only valid in the absence of finite size ef- 
fects and suggested that finite size could result in large 
values of z&3 

Finite size effects, which include a finite size-system 
and the inherent magnetic penetration depth \±£3 in su- 
perconductors, have long been known to complicate the 
verification of KTB behavior because they make free vat 
tices possible for any temperature range. Kadin et a/.EJ 
and Lee and Garlands were among the first to study 



this, experimentally while Repaci et alS3 and Herbert et 
al.l3 have recently extended this work. Theoretically, 
Simkin and KosterlitztH have investigated numerically 
and analytically, the finite size effects in these systems, 
while ArtemovLa has examined the effect of "Pearl" vor- 
tices on the magnetization in superconducting thin films. 
And more recently, Pierson and Validly have completed 
a renormalization group study of the two-dimensional 
Yukawa Gas, which incorporates the inherent finite size 
Aj_ into the two-dimensional Coulomb Gas (2DGG). 

Finite size (FS) scaling, which incorporates the effect 
of a finite size into the FFH scaling, is an established 
method in tho-study of numerical data in KTB behavior. 
Lee and TeitelEl were among the first to use this approach 
to show that the dynamic critical exponent z has a value 2 
for the two-dimensional Lattice Coulomb Gas (2DLCG). 
Medvedyeva et alS3 have also used finite size scaling to 
reinforce the result z = 2 for the 2D resistively-shunted- 
junction model. (Numerical studies in fact have provided 
many insights into KTB behavior. A review of those 
results is beyond the scope of this paper but we re£er.-the 
reader to many excellent papers in the literature.tjEj) 

Like the conventional approach and FFH scaling, FS 
scaling has its strengths and weaknesses. For example, 
because FS scaling requires the ability to carefully tune 
the system size, it is not conducive to studying experi- 
mental data where the finite length scales are more illu- 
sive. It is therefore important to keep in mind that all 
approaches provide a great deal of information and often 
complement one another. 

The correlation length is an important quantity in FFH 
scaling but questions about its form below T c remain. In 
the original paper of Kosterlitz,0 the correlation length 
was taken to be infinite below the transition temperature 
owing to the divergence of the-.susceptibUity there. In 
subsequent papers by Friesen,E3 Pierson,E3 and Simkin 
and KosterlitzjUl the vortex correlation length is taken to 
be finite for this temperature regime (but still diverging 
at T c from below) and interpreted roughly as the size of 
the largest vortex pairs. Other authorst3 maintain that 
the correlation length is infinite for T <T C . 

In this paper, we do a numerical study of the 2DLCG 
using Monte Carlo simulations to understand the differ- 
ence in the value of z obtained from finite size scaling and 
FFH scaling. We vary the system size to study the effect 
of finite sizes and verify that z = 2 using FFH scaling, 
thereby making FFH scaling consistent with FS scaling. 
We also analyze the transport characteristics and find ev- 
idence that the vortex correlation length is finite below, 
the phase transition, in contrast to the claims of others.Ej 
The numerical approach that we use is that of Lee and 
TeitelEfi extended to larger system sizes [L = 200). 

The paper is organized as follows: Section II describes 
the model that we use along with the numerical and scal- 
ing methods. In Section III, we present our results and 
analysis of the transport characteristics. Implications 
and a summary of our work is presented in Section IV. 



II. MODEL AND METHODS 

A. 2D Lattice Coulomb Gas with an applied field 

To study the KTB critical behavior and its dynamic 
critical exponent, we use the 2DLCG model, whose re- 
lationship to the KTB transition and mapping to the 
2DXY model are well knownHffia This model is described 
in Refs. 
here. 



20 and E2| and so will be only briefly described 



The Hamiltonian for the 2DLCG is given by 



H = -x^2,qiV{r i:j )q h 



(1) 



where the sum is over all pairs of sites in the lattice, q^ 
is the (integer) charge at site i, r^ = 



Tj, and r. 



is the location of the ith charge. V, which solves the 
discrete form of Poisson's equation in two dimensions, is 
the Coulomb potential (and is therefore logarithmic at 
large distances). Using Fourier tranforms, V in Eq. (fil) 
may be replaced by 



V(r) 



S?™ 



ik-r 



where, for the square lattice considered here, 



2 — cos(k ■ x) — cos(k ■ y) 



The wave vectors are defined in the usual way: 
2irmi , 27TTO2 „ 



L 



L 



(2) 



(3) 



(4) 



where mi, m^ = 0, 1, 2, . . . , L — 1. In addition, since we 
wish to apply our results to 2D superconductors in zero- 
magnetic field, we use the constraint that the total charge 
in the system must be zero, which deals with the infrared 
divergence as k — ► 0. 

In the presence of an electric field, E, an electric po- 
tential energy teem of the form J^ q^i ■ E is added to 
the Hamiltonian.EJ The new term results in a net move- 
ment of charge along the direction of the field E that we 
measure in terms the charge current density, J, using the 
approach of Ref. |2^. 

The applied field E and the resulting vortex current 
density J in the 2DLCG model are analogous to the cur- 
rent I and voltage V in a superconducting thin film, but 
in an inverted way due to the mapping between the two 
systems. In the superconducting thin film, the current 
is injected into the sample, which causes vortex move- 
ment in a direction perpendicular to the current via the 
magnus force. This results in resistance and therefore a 
voltage V. For the 2DLCG, it is the voltage that is ap- 
plied and that results in the vortex "charge" current in a 
direction parallel to ther-electric field. Following the no- 
tation of Lee and Tcitcl,E2l the following correspondences 



between the two systems apply: E < — ► I and J < — ► V. 
In this paper, we will examine the E-J curves of the 
2DLCG, which correspond to V-I characteristics in ac- 
tual superconducting thin films. Unless otherwise noted, 
the notation that we use will be that of the 2DLCG. 
Along these lines, the E-J curves will be plotted as E 
versus J in our figures, so that they mirror the analo- 
gous curves in the 2D superconcuctors that are plotted 
as V versus I. We will therefore use the notation V-I for 
the current-voltage curves instead of the more standard 
I-V. 



where a smaller number of MCS was used due to very 
long run times. At high electric field, the current is large 
and the shorter runs provide sufficient accuracy. At each 
E and T, at least five runs were performed. The error 
bars on the points in Figure EL Figure 0, and Figure H 
are the standard deviations of the average current from 
these runs. In some cases at low E, many additional runs 
were performed in an attempt to reduce the error bars. 



C. Fisher-Fisher-Huse and Finite-Size Scaling 
Methods 



B. Monte Carlo Simulations 

The Metropolis^ Monte Carlo technique is used to cal- 
culate the J-E data based on the Hamiltonian discussion 
in Section IIA. Briefly, the technique, as it applies to the 
2DLCG model with a square (N = L X L) lattice and 
integer charges, is as follows: (See Ref. g0| for more de- 
tails.) 

1. Create a trial system from the original system 
by choosing a site, i, and one of its four neighbors, 
j, at random. Add one unit of charge to site i and 
subtract one unit from site j so that the net charge 
in the system is always zero. 

2. Compute AH = H trial - H original . Since only 
two sites have been altered, most terms in H tria i and 
Horiginal are the same, simplifying the calculation of 
AH. 

3. If AH < or if e - AH / T is greater than a random 
number uniformly distributed on the interval [0,1), 
accept the change and replace the the original sys- 
tem with the trial system for subsequent trials. Oth- 
erwise, leave the system unchanged. 

When evaluating the Hamiltonian for a trial system, 
the minimum image convention and toroidal (periodic) 
boundary conditions are used.M In addition, a "stability" 
term of the form Yli(li~ 1i) ^ s added to the Hamiltonian 
to minimize the probability of creating multiply charged 
sites in the system. E3 

For many combinations of E and T, particularly when 
the charge density is small, a low rate of acceptance of 
trials is expected. When this occurs, it is convenient 
to calculate AH without summing over the array. The 
methO|d-| employed here, used and described by Lee and 
TeiteljH3 is to store information on each site in the system 
in the form of the potential, which is updated for each 
site after each accepted trial. This technique can improve 
run times on the order of the reciprocal of the acceptance 
ratio, often factors of O(10 3 ) or more. 

A sequence of N trials is called one Monte Carlo Step 
(MCS). All of our computer runs consist of an initial 
10, 000 MCS to equilibrate the system during which no 
data are saved, followed by at least 500, 000 MCS, for 
which the average current is computed. The only ex- 
ceptions were some large systems (L = 200) at E > 0.6 



Given the difficulties of the conventional approach to 
analyzing V-I characteristics discussed in Section |l it is 
important to use other methods including the powerful 
scaling methods. In this paper, we will use Fisher-Fisher- 
Huse scaling and finite size scaling. 

In the FFH theoryB for 2D superconductors, the E-J 
(in the notation of the 2DLCG) curves should scale as 



J = EC'X±(E£/T), 



(5) 



where X+(—) ( x ) i s the scaling function for temperatures 
above (below) Tktb- The two important asymptotic 
behaviors of x( x ) are lininoX+t 1 ) = const. (Ohmic 
limit), and linix^oo x± ( x ) °^ Smi-tepifical isotherm). As 
established in previous papers JatJJ'liil it is advantageous 
to rewrite Eq. @ as 




= £± (££/T) 



(6) 



1 I 7 

(where £±(x) = x/x± (x)) because only the x-scale is 
stretched in Eq. (||). In Eq. (||), both the a;-scale and 
y-scale are stretched making it harder to judge a collapse 
of the scaled data. 

The FFH scaling has three fitting parameters z, T c , 
and b, which can be determined in the three scaling steps. 
The first, which yields estimates of z and T c , is to deter- 
mine the critical isotherm. In the absence of finite size 
effects, the critical isotherm is straight on a log-log plot 
since J ex E z+1 for all E at T — T c . The second step is 
to estimate b in the exponent of the correlation length, 
which, for experimental data, one can do by examining 
the resistance since R(T) ex £~ z . EoEJiumerical data, one 
assumes £ ex exp[C/|T — Tcl 1 / 2 ]!^!^! which is generally 
expected to hold for T > T c . The final step is to achieve 
a scaling collapse for Eq. (0), adjusting the three param- 
eters within the error bars established in steps 1 and 2. 
One should also bear in mind that all these steps can be 
repeated to optimize the scaling collapse. 

In the finite size scaling form, valid for £ > L, the 
system size L is the dominant length scale_and therefore 
substitutes for £(T) in Eq. (||) to arrive at£3 



JL Z /E = X ±(EL). 



(7) 



Strictly, this equation works only for T = T c , but is also 
valid for T ~ T c since leading order corrections are small: 
0(L/£). In this powerful approach, there is only one 
fitting parameter, z, making its results more definitive. 

To achieve the best scaling collapse for the FFH and 
FS scaling, the principle method that we use is to sort 
the data according to the value of the scaling variable 
x (whether it be E£/T or EL) and then to minimize 
the difference of the logarithm of the values of two. ad- 
jacent, scaled points using a least squares method!] As 
a second means of optimization, we tried approximat- 
ing the T > T c curve to a function and then minimizing 
the difference between the scaled points and the curve. 
However, because of the large number of fitting param- 
eters, definitive results could not be achieved with this 
approach. 

For various reasons, not all of the data are used in the 
optimization. At large values of the applied electric field 
(E ~ 0.8), the isotherms become ohmic due to a satu- 
ration of the vortex density. (See discussion in Section 
[IIA.) Because this behavior is not due to critical behav- 
ior, the data should not scale (collapse) with the other 
data and so is not included in our optimization. How- 
ever, we do plot those points in the scaling plots (Figures 
even though those points do not collapse. (In 
10| and [by, the analogous, large / data was not 
plotted with the scaled data.) In addition, because the 
dynamic critical exponent z has a temperature depen- 
dence below T c due to the line of fixed points in the KTB 
critical behavior, the T <T C data are not used in the op- 
timization of the FFH scaling collapse. This is because 
the temperature dependence of z is not the same as the 
scaling variable x — E£/T and so the FFH .scaling is not 
expected to be as effective in this regime.Elc3 Nonethe- 
less, we do scale the T <T C data assuming a symmetric 
correlation length. Finally, data are also excluded from 
the optimization if the error bars are greater than 50%. 



III. RESULTS 

The principal objectives of our work are to study the 
value of the dynamic critical exponent z for the 2DCG 
in the context of Fisher-Fishcr-Husc scaling and the in- 
fluence of finite size effects, and to study the behavior of 
the correlation length £_(T) below the phase transition. 
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FIG. 1. Charge current density J(E,T) versus E for fixed 
system size L = 24 for temperatures ranging from 0.09 to 0.25 
by increments of 0.01 (from bottom to top). 



A. Dynamic Critical Exponent z 

The two principal system sizes used in our simulations 
were L = 24 and L = 100. The J-E isotherms for the 
L = 24 system are shown in Figure |j. As one can see the 
isothermal J-E cpsves mimic closely the V-I of thin-film 
superconductors.Ej Reflecting the absence of vortex pairs, 
the higher temperature curves are primarily Ohmic (i.e., 
have slope 1 on a log- log plot). The lowest temperature 
curves are primarily non-Ohmic, indicating the presence 
of vortex pairs. The intermediate-temperature curves, 
on the other hand, are Ohmic for small E, become non- 
Ohmic at larger values of E, and abruptly turn Ohmic 
for the largest values of E. This of course reflects the 
fact that with each different value of E, one is probingj-a. 
different length scale, which varies as the inverse of £\l3 
(See discussion in last paragraph of Section IIIB.) Conse- 
quently, at small E, long length scales are being probed 
where free vortices exist and Ohmic behavior is expected 
for T > T c (or when Aj^ is finite). For large values of E, 
smaller length scales are probed where vortices are still 
bound in pairs and result in non-Ohmic behavior. The 
other salient feature of this graph is that the error bars 
become large at the small- E end of the isotherms due to 
small numbers of vortex pairs. 
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10 10000 
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24 data of Figure hi scaled according 
i) using z = 6.01, T c = 0.127, and 



The temperature of the phase transition in the L = 24 
data is not obvious in Figure [y. Using the FFH dynamic 
scaling criteria for the critical isotherm, which says that 
J oc E z+1 for all values of E at T = T c , results in T c ~ 
0.14 and z ~ 6. The subsequent FFH scaling for the 
T > T c data yields an excellent collapse for z = 6.01, 
T c = 0.127, and C = 0.7475, as seen in Figure |. (As 
noted in Section ||, the crossover to Ohmic behavior at 
large E data seen Figure h is not due to critical behavior 
and so should not be used to judge the collapse. It is this 
large- E data that does not collapse with the other data 
in Figure 0.) 



t = 0.170 



t = 0.200 



t = 0.220 t = 0.250 

FIG. 3. Sample vortex configurations for E — 0.3 at 

four temperatures. A white box is a vortex, a black box an 

anti-vortex, and a grey box an empty lattice site. 



The values of the parameters in this collapse, opti- 
mized using our method described in Section [Il|, are con- 
sistent with the lawre-.values of z and low values of T c 
found previously.clil3 , Lj However, an examination of thg 
inverse dielectric constantE2l e _1 and the vortex densityEl 
reveals no critical behavior at that temperature with or 
without a current. An inspection of the lattice config- 
urations also reveals no clues of a phase transition at 
T ~ 0.13. In Figure |l characteristic vortex configu- 
rations in the lattice are show for various temperatures 
for a relatively large value of the electric field, E = 0.3. 
Even at the lowest temperature shown here (T = 0.17), 
there are few vortices and they are clearly bound to one 
other. (One can also see in that figure that the vor- 
tices seem to become unbound above the expected T c of 
0.218.) Hence, the apparent FFH scaling collapse at the 
large value of z and the small value of T c does not seem 
to correspond with a real phase transition. We now in- 
vestigate if the FFH scaling collapse for these values of 
T c and z could be due to finite-size effects, as recently 
suggested by Medvedyeva et al£3 
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FIG. 4. The J-E curves for T = 0.210 for a variety of 
system sizes. [Inset: The same set of data scaled with finite 
size theory (Eq. m) for z — 2.15.] 

To further explore the value of z and the transition 
temperature, we employed finite-size scaling methods for 
T = 0.210, ^temperature just below the expected value 
T c = 0.218.E3L3 The J-E curves for values of L ranging 
from 12 to 200 are shown in Figure H. As one can see, 
the curves become straighter (on a log-log scale) as ex- 
pected near the critical isotherm when the system size is 
increased, indicating that the low-_B Ohmic behavior for 
the L = 24 data is attributable to finite size effects and 
that the true critical isotherm occurs at a larger tem- 
perature. The finite size scaling, shown in the inset to 
Figure [|, reinforces this idea. With z being the only 
fitting parameter, the optimal value is, 2. 15, reinforcing 
the original results of Lee and Teitclc2l (and the subse- 
quent work of others) that z — 2 and that T c occurs at a 
much larger value than obtained from the FFH scaling. 



The fact that the value of z here is slightly larger than 
2 can be attributed to the fact that the temperature is 
just below that of the expected value of the transition 
temperature and that z increases in value for T < T c . 
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FIG. 5. Charge current density J(E, T) versus E for fixed 
system size L = 100 for the temperatures 0.09, 0.12, 0.15, 
0.16, 0.17, 0.18, 0.19, 0.195, 0.20, 0.205, 0.21, 0.211, 0.212, 
0.213, 0.214, 0.215, 0.22, 0.23, 0.24, 0.25, 0.28 (from right to 
left). The density of curves is larger near T c . 

The value of z obtained in the finite size scaling re- 
inforces the idea that the large value of z found for the 
L = 24 data in the FFH scaling is due to finite size effects. 
We have therefore calculated J-E curves for L = 100. 
(See Figure g|.) The J-E curves are similar to those 
of L = 24 both in terms of qualitative behavior and in 
terms of the error bars. The primary difference is that 
the intermediate-temperature curves are straighter than 
those of L — 24 and the critical isotherm now occurs at a 
temperature T ~ 0.210 for z ~ 2. (Of course, if we could 
extend the J-E data to smaller J, the J-E curves would 
not be straight for T < T c due to finite size effects.) Us- 
ing these values as a guide, the best collapse for this data 
with the FFH scaling, shown in Figure |6|, was obtained 
using z = 1.90, T c = 0.212, and C = 0.65. As one can 
see, the collapse is quite good for the T > T c data. (Note 
that four points for T = 0.22 were excluded from this 
figure due to the fact that the error bars were greater 
than 50%.) Even the T < T c data scale well despite 
the fact that z in this regime has a temperature depen- 
dence different than that of the scaling parameter x. The 
T = 0.090 and T — 0.12 are the exception to the collapse 
but is most likely due to the fact that that data is well 
outside the critical region. The fact that the rest of the 
T <T C data scales reasonably well may be a coincidence 
or evidence for a finite correlation length. (We present 
evidence for a finite correlation length for T < T c in the 
next subsection.) The value of the T c achieved through 
our methods is slightly smaller than the expected value, 



which we attribute to finite size effects. We add that an 
acceptable collapse is also achieved for T c — 0.218 and 

z = 2. 




1000 100000 
E^/T (arb. units) 

FIG. 6. The L = 100 data of Figure ^ scaled according to 
FFH scaling (Eq. % using z = 1.90, T c = 0.212, and C = 0.65. 
The T = 0.09 and T = 0.12 isotherms do not collapse with 
the other curves, which we attribute to those curves being far 
from the critical region. 

We believe that this work is a definitive verification 
of z = 2.0 ± 0.2 found using FFH scaling and therefore 
puts FFH scaling in accord with finite size scaling. We 
conclude that the larger valuc-of-,z (5.6) found for the 
L = 24, data and elsewhereoEIlO is due to finite size 
effects.EHl (See discussion in Section IV.) This work also 
vindicates the use of FFH scaling for systems free pf finite 
size effects. Lee and TeitelHJ and Harris et als£ have 
previously showed a, z — 2 collapse using FFH scaling for 
their numerical data (L = 24) and JJA data, respectively. 
However, we do not believe that that work should not be 
taken as evidence for z — 2 because the scaling collapse 
was not convincing, they did not vary z to optimize the 
collapse, and it was later shown that the data collapse 
better for much larger values of z.Qil 

Finally, we note that we calculated the values of the 
J-E exponent ajn order to study the predictions of 
Minnhagen et alEB about its temperature dependence. 
In accord with the work of others J3EJ our work indicates 
that the prediction of Ref . g5| for the temperature depen- 
dence of a (usually denoted cvpm) is in better agreement 
with the data than the original prediction^ for the tem- 
perature dependence of a (sometimes denoted ocahns-) 
This agreement further verifies our L = 100 J-E data. 



B. Vortex Correlation Length for T < T c 

Our method for studying the vortex correlation length 
is to investigate the value of the J-E exponent a as a 
function of E by looking at the slope of our J-E curves. 
More specifically, we study the derivative dlog J/d\ogE, 
which is equal to a when J oc E a . A similar examination 



by Medvedyeva et alE-z (see their Figure 2) for system 
sizes of L = 8, using the analogous V-I curves in the 2D 
RS J model, revealed that the slope started at a low value 
for small J, peaked, and then descended to a value near 1 
at large values of I. For each of the isotherms with tem- 
peratures below the critical temperature, they-Jfound the 
peak to occur at roughly the same value of /Jia Because 
the scaling variable is either I£ when ^Ctor IL when 
L < £, the fact that the peak in d log V/d log I does not 
move is an indication that finite size events dominate this 
temperature regime (§*id the scaling parameter is EL). 
Above T c , the authorsI3 found the peak to move quickly 
to larger values of I, indicating that I£ is the princi- 
pal scaling variable since £ is decreasing quickly in this 
regime. 




FIG. 7. The L = 24 J-E data plotted as dlog J/dlogE 
versus E for all of the temperatures shown in Figure 111 The 
top curve corresponds to T — 0.09 and the bottom curve to 
T = 0.25. 

For our L — 24 data, d log J/d log E behaves differently 
below T c than does the analogous quantity of Ref. O. As 
one can see in Figure R, the locations of the peaks as a 
function of E move to smaller values as one increases the 
temperature from 0.09 to about 0.14. In the temperature 
range 0.15-0.22, the peaks develop into plateaus. 

The peaks and plateaus indicate that there are two 
competing length scales. On the large length scale 
side (small E) is the temperature-dependent correlation 
length £_ (T) while at small length scales (large E) is the 
lattice size, which enters because of a saturation of vortex 
density. A plateau signifies that £ is much larger than the 
saturation length scale while a peak indicates that £ is 
approaching the lattice size. The value for the left hand 
size of the plateau is determined by £ or L, whichever 
is smaller. In Figure R, the point at which £ becomes 
larger than L occurs at about T = 0.19, above which 
the plateau begins at E ~ 0.13 on the low E side. That 
the value of E for which d log J/d log E peaks or begins to 
plateau depends on temperature for a broad temperature 
range below T c is an indication that the scaling variable 
for T < 0.19 < T c is E£_ and that, therefore, £_(T) is 



finite below T c . 

The lattice spacing affects the J-E curves for E > 0.7, 
as can be seen in Figure m, and determines a maximum 
value of E at which dlog J/d log E peaks. Therefore, the 
discrete lattice size can also result in the position of the 
peak in dlog J/d log E not moving. As a consequence, 
that the peaks would be located all at the same approxi- 
mate value of E, as found in Ref. [lj, should not be taken 
as evidence strictly for finite size effects. 

The values of E that correspond to the influence of 
the lattice spacing and the lattice size can be predicted 
by simple arguments. (See for example Ref. ||.) The 
approximate energy for a vortex pair is (logi? — ER)/e, 
which increases for small R and then decreases for large 
R, thereby peaking at R = 1/Ex5 (This equation ex- 
plains why one probes different length scales by changing 
the value of E as discussed in Section III A.) Since the lat- 
tice size corresponds to a value 1, one expects its effects 
to manifest themselves for E ~ 1, which is what we ob- 
serve in Figure 0. Similarly, one expects finite size effects 
to manifest themselves for E ~ 1/L ~ 0.04. Indeed, in 
Figure R, all of the isotherms are ohmic for E < 0.5, in 
agreement with our prediction. 



IV. DISCUSSION AND SUMMARY 

Using system sizes up to L = 200 in Monte Carlo sim- 
ulations of the 2DCG, we conclude z = 2.0 ± 0.2 us- 
ing Fisher-Fisher-Huse scaling and finite size scaling, re- 
inforcing the results of others that z — 2. (See, e.g., 
Refs. |l2|)P0|,pi| .) We have also presented evidence that 
£_, the correlation length below T c , is finite. 

Previous analysis of experimental data on 2D super- 
conductors, JJA's, and superfluids, androf-jpperical data 
on the 2DCG by two of these authorseNo using FFH 
scaling indicated z ~ 5.6. As stated here and in Ref. Ill, 
it appears that this large value is due to finite size effects. 
Because a majority of the V-I data from the literature 
was examined in Ref. 0, it appears that nearly all exper- 
imental V-I measurements on 2D superconductors and 
JJA's are dominated by finite size effects. Why the value 
Z ~ 6 should result from finite size effectsEj remains a 
mystery, and should be studied further. 

Because these results indicate that finite size effects 
are ubiquitous in 2D superconductors and 2D JJA's, the 
V-I and R data from the literature that were used to 
verify KTB behavior in the conventional manner needs 
to be re-examined to account for the finite size effects. 
In the presence of strong finite size effects like those in 
these samples, two aspects of the conventional method 
fail. First, the small / behavior of the V-I isotherms 
would not reveal any quick changes in a near T c because 
finite size effects dominate at small /. In other words, as 
is well known, finite size effects wash out the expected 
jump in a. Secondly and more importantly, the resis- 
tance formulas are no longer valid in the presence of fi- 



nite size effects. That previous authors could fit their 
resistance data to the more commonly used Minnhagen 
formula is most likely due to the number of fitting param- 
eters. Finally, we stress that the most effective means for 
studying the temperature dependence of a should begin 
by extracting a(T) from the V-I curves in 2D supercon- 
ductors and JJA's for the same value of /. 

Clearly, in order for this data to be re-interpreted, a 
better understanding of finite size effects is needed since 
previous checks for finite size effects in the data of Ref . [15| 
were negative. For example, those authorscJ found that 
the value of Aj_(0) derived from a fit to Aj_(T) was four 
times larger than, the expected value. Furthermore, for 
that same data,ll3 where finite size effects are now known 
to dominate, the resistance was found to fit better to the 
Kosterlitz resistance formula than a resistance formula 
derived for finite-size dominated Ohmic resistances Fi- 
nally, in Ref. |[ it was found that the length scale at 
which the crossover from finite-size-induced Ohmic be- 
havior to non-Ohmic behavior was expected to occur for 
the data of Ref. [15] did not correspond with the observed 
value. 
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